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Electrons/atoms can flow without dissipation at low temperature in superconductors/superfluids. 
The phenomenon known as superconductivity/superfluidity is one of the most important dis- 
coveries of modern physics, and is not only fundamentally important, but also essential for 
many real applications An interesting question is: can we have a superconductor for heat 
current, in which heat/energy can flow without dissipation? Here we show that heat super- 
conductivity is indeed possible. We will show how the possibility of the heat superconduc- 
tivity emerges in theory, and how the heat superconductor can be constructed using recently 
proposed time crystals 0^1 The underlying simple physics is also illustrated (see discussions 
after Eq. @). If the possibility could be realised, it would not be difficult to speculate various 
potential applications, from heat tele-transportation to cooling of information devices. 

The usual superconductivity emerges as a result of spontaneous breaking of electromagnetic 
gauge symmetry 4 . The gauge symmetry dictates that all electron systems are invariant under 
changes of phase pre-factors of wave amplitudes of constituent electrons. The conservation of 
electron numbers is a result of the symmetry, according to Noether's theorem. A superconductor 
breaks the gauge symmetry: its ground state is not a phase symmetric and fixed-particle-number 
state, as expected from the gauge symmetry. Instead, a macroscopic superconductor has definite 
phase for its constituent particles (Cooper pairs). When coupling two superconductors with differ- 
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ent phases, one finds a super-current between them, known as Josephson effect 1 . More generally, 
a slow spatial variation of the phase in continuous space will drive a super-current. 

The above discussion provides useful hints on how the heat superconductivity could emerge. 
Like the electron numbers, energy is a conserved quantity. Actually, it is the only quantity that is 
conserved for all imaginable physical systems. The symmetry associating with the conservation of 
energy is nothing but the time translational symmetry. An analogy to the usual superconductivity 
immediately suggests that the heat superconductivity could emerge as a result of spontaneous 
breaking of the time-translational symmetry, which is equivalent to spontaneous emergence of 
local clocks^. One would expect a spatial variation of local clocks could drive a heat super-current, 
i.e., the Josephson effect for heat. 

The time-crystal proposed by Wilczek et al. provides an instance of spontaneous time trans- 
lational symmetry breaking 03. A couple of model systems to realise the time crystal had been 
proposed, including soliton phase of a collection of charge particles with attractive interaction 
confined in a ring ^, and Wigner crystal phase of a group of ultracold ions confined in a ring- 
shaped trapping potential 3 . In both cases, the system undergoes spontaneous symmetry breaking 
of spatial translational symmetry, and forms a phase of inhomogeneous density distribution. As a 
result, there emerges a new parameter Oq, which specifies center of the soliton or the origin of the 
Wigner crystal. In the macroscopic limit of the particle number N — > oo, there is emergent orthog- 
onality between different ^-states, and for all local observations, the system could be considered 
to be in a state with definite 6q ^. 
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A magnetic field transforms the system to a time crystal, which can be visualised as a persis- 
tent rotator with a finite angular velocity in its ground state 

u = (l + a)j, (1) 

where a = qcf)/h with <p being the total magnetic flux threading through the ring, q is charge of 
ion/particle, I is momentum of inertia of an ion/particle relative to the center of the ring, and I is an 
integer (or a half-integer for an odd number of fermions) that minimises \l + a\. The corresponding 
ground state wave function can be related to the zero magnetic field 6> -state 6 ) by^ 

= e~ il ^ e ^({e i },9 -u J ot), (2) 

where {#;} denotes the set of coordinates of particles on the ring. When a is not an integer or a 
half-integer^, the ground-state many-body wave-function becomes time-dependent, and the system 
becomes a persistently rotating object. Here, 9 acquires a new meaning: the phase angle of the 
persistent rotation. Effectively, it defines the local clock. 

A heat superconductor can be constructed from a group of time crystals by stacking them 
concentrically to form a three-dimensional rod, as shown in Fig.[TJ At sufficiently low temperature, 
the inter-ring coupling due to the density inhomogeneity will stabilise a configuration of relative 
phase angles of the rings. It also gives rise to a finite torsional stiffness. The elastic energy induced 
by a twist of the rod can be written as, 

E-E =^j dx[d x 9 (x)] 2 • (3) 

in the continuum limit, where E is the ground state energy, k is the torsional stiffness, and we 
assume that the rod is extended along the x-direction. 
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To determine heat current along the rod, we introduce a gravitational vector potential A E {x), 
which couples to the heat current in a way just like that of the usual vector potential to the electric 
current. We apply the local gravitational gauge symmetry (see Methods): the system should be 
invariant under joint transformations (x,t) — >■ (x,t + Xe{x)) and A E (x) — > A E (x) + 8 x xe(x) 
to the first order of a slowly varying function Xe(x). For the ground state wave function Eq. ^ 
of time crystal, the local time translation is equivalent to a translation of the phase angle 9q{x) — > 
6q(x) — cuoXe(x). The symmetry dictates the minimal coupling form to the gravitational vector 
potential: 

E-E = ^Jdx [dj (x) + u A E (x)] 2 . (4) 



The energy current along the rod can then be determined: 

= -Kuj d x 9 (x) . (5) 
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We see that a heat (energy) super-current emerges in a twisted rod in the presence of the persistent 
rotation uq. The result is consistent to our general theoretical consideration, as a spatial variation 
of ^o(^) implies a variation of the local clocks. 

Equation Q is the central result of this paper. The underlying simple physics is illustrated in 
Fig. [2] We consider a heat superconducting rod in contact with environment through two viscosity 
fluid reservoirs. Had the centre rod been the usual solid, the viscosity fluid would exert friction 
torque that will eventually put it to rest. However, a rod made of the time crystals rotates persis- 
tently at the angular velocity u , and cannot be slowed down because it is already in the ground 
state. 



We consider the situation when the environment exerts torques Tl(r) at the two ends of the 
rod. In the presence of the persistent rotation, the torques will do works Wl(r) to the rod^ 

dW fa R) = u oTl(r) ■ (6) 

The required energy must be absorbed from or released to the environment. At the stationary state 
that the phase angle configuration of the rod is fixed, there is no net energy exchange between the 
rod and the environment, and we must have: 

(Tr) = -(Tl)=T. (7) 

As a result, the energy will be absorbed from the environment at the one end, and released to the 
environment at the other end. Effectively, a heat current uqT passes through the rod from one end 
to the other. At the same time, an average torque T will induce a twist: 

d x e {x) = -T/K. (8) 

Combining Eq. ([6]) and ([8]), we obtain the heat super-current formula Eq. ([5]). 

The emergence of the torques at the two ends of the rod is resulted from the microscopic 
process of energy exchanges with the environment. At the finite temperature, the two ends of the 
rod continuously receive random torsional kicks from the environment. Each random kick pushes 
the system out of its ground state temporarily, accompanying with energy exchanges. The system 
reaches equilibrium when energy absorptions and emissions are balanced. For a time crystal rod, 
the process gives rise to temperature dependence of equilibrium expectation value of the persistent 
angular velocity u)q(T) 3 . When the system deviates from the equilibrium, there will be net energy 
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influx or outflow, corresponding to an average torque: 



where J E is the net energy flux. 



We now consider an infinitesimal temperature difference 5T such that the left (right) reser- 
voirs in Fig. [2] has temperature T L ^ — T ± 5T/2. For a symmetric structure, the rod must be 
at the temperature T, maintaining a persistent angular velocity ujq(T). It is important to note that 
the whole rod must be at the same temperature so that the different parts of the rod could rotate 
synchronously. The temperature differences between the rod and the environment at the two ends 
will drive net energy fluxes, 

J L E (R) = ±4 n y > do) 

where we introduce an interface heat conductance for heat conduction between the normal 
and superconducting heat mediums. The corresponding torques determined by Eq. ([9]) twist the 
rod, inducing heat super-current in the rod. The calculation of k^ n depends on details of mi- 
croscopic model, in analogy to the calculation of the tunnelling conductance of a normal metal- 
superconductor junction in the usual superconductivity H To transport heat, one does need a tem- 
perature difference to overcome the contact heat resistance 2/k^ n of the two interfaces. On the 
other hand, the heat flow generated by a given temperature difference does not depend on the length 
of the central heat superconducting rod, meaning that one can transmit heat for arbitrary distance 
without incurring transmission loss in the rod. 

There is also AC Josephson effect counterpart in heat superconduction ^ This can be 



achieved by cutting the centre rod of Fig. [2] into two weakly coupled disjointed pieces. The two 
pieces will have different persistent angular velocities ojq(Tl) and Uo(Tr) when the temperatures 
of the two reservoirs are different. The resulting heat super-current will oscillate at frequency 
\u (T L ) — u (T R ) \ or one of its harmonics. 

In summary, we show that the heat superconductivity naturally emerges in a system spon- 
taneously breaking the time translational symmetry. Although our discussion is based upon the 
existing models of time crystal and has a mechanical intepretation, the essential physics is com- 
pletely general, i.e., a spatial variation of local clocks will drive a heat super-current. The pos- 
sibility to get a real heat superconductor in laboratory is certainly rested on the realisation of the 
time crystal, which is still challenging both theoretically and experimentally. On the other hand, 
one should not ignore the possibility of finding the heat superconductivity in natural compounds, 
as the usual superconductivity was first discovered in mercury. Nevertheless, it is yet to see any 
fundamental principle prohibiting the spontaneous breaking of time translational symmetry from 
happening, and a connection between it and a new matter of vast application potential (i.e., heat 
superconductor) would be a great incitation for future pursuits. 

Methods 

Luttinger introduces a gravitational scalar potential ip(r) which perturbs the hamiltonian^: 

H-> I dr [1 + ip(r)] h(r) , (11) 
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where h(r) is the local hamiltonian density, and H = J drh(r). One can further introduce gravi- 
tational vector potential A E {r) which is coupled to the energy current operator j E (r): 



H -> J dr [h(r) - j E {r) ■ A E (r) 



(12) 



The set (ip, Ae) provides heat transport counterparts of electromagnetic scalar and vector poten- 
tials HE. 

The energy current operator j E (r) is defined by the energy continuity-conservation equa- 
tion^ 

° hlr) ' \h(r),H]=-V-j E (r). (13) 



dt \h 

The equation only defines the energy current up to a curl. To eliminate the uncertainty, we further 
impose the scaling law^^ 

j E (r) -»- [l + ^(r)] 2 j E (r) (14) 



when the hamiltonian is scaled by the gravitational scalar potential as Eq. (11). 



We can introduce the local gravitational gauge transformation: 



Uq = exp 



i/h) / drx E (r)h(r) 



(15) 



where XE( r ) is a slowly varying function. It is easy to see that the transformation is just a local 
time translation (r, t) — >• (r, t + Xe{i*)) if one ignores the gradient corrections of xe(cc). 

Under the local gravitational gauge transformation, the local hamiltonian density is trans- 
formed, to the first order of Xe{t): 



h(r), / drxE{r)h(r) 



h(r) + VxE-j E (r) + V' 



Xe3e\t i 



■ (16) 



The commutator can be derived from the scaling law Eq. ( 14) 10 . 



Using Eq. ( |T6| ), it is easy to verify the local gravitational gauge symmetry: to the first order 
of XE^f), the total hamiltonian H is invariant under the joint transformations Uq and Ae(v) — > 
A E (r) + V XE (r). 
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Figure 2 
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Figure 1 (a) The structure of a heat superconducting rod. A series of soliton rings/Wigner 
crystal rings are arranged concentrically along the direction perpendicular to the rings. 
The rings are stabilised to a configuration of relative phase angles by the inter-ring cou- 
plings. All rings rotate synchronously and persistently at angular velocity u . A conducting 
filament is threaded through the central axis of the rod and in contact with external elec- 
trodes, providing free charges to compensate the charges trapped on the rings. Here, the 
particle density profiles of soliton rings are shown and colour-coded. The spacings be- 
tween rings are exaggerated, (b) A twist of the rod induces a heat super-current, indicated 
by the magenta arrow. 



Figure 2 A heat superconductor in contact with environment. A structure of Fig. [T|is as- 
sembled and represented by a rod in the middle. The rod has a persistent angular velocity 
u . The two ends of the rod submerge into two reservoirs of viscosity fluid (represented 
by two cylinders at the two ends). The directions of friction torques exerted by the fluid 
reservoirs to the rod are shown by the blue arrows. In the presence of both the torques 
and the persistent rotation, there will be energy exchange between the rod and the two 
fluid reservoirs. The corresponding directions of the energy fluxes are shown by magenta 
arrows. At the stationary state, the friction torques at the two ends must have the same 
magnitude and the opposite directions. 
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